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ABSTRACT 

21-cm emission from neutral hydrogen during and before the epoch of cosmic reioni- 
sation is gravitationally lensed by material at all lower redshifts. Low- frequency radio 
observations of this emission can be used to reconstruct the projected mass distribu- 
tion of foreground material, both light and dark. We compare the potential imaging 
capabilities of such 21-cm lensing with those of future galaxy lensing surveys. We 
use the Millennium Simulation to simulate large-area maps of the lensing convergence 
with the noise, resolution and redshift-weighting achievable with a variety of idealised 
observation programmes. We find that the signal-to- noise of 21-cm lens maps can far 
exceed that of any map made using galaxy lensing. If the irreducible noise limit can 
be reached with a sufficiently large radio telescope, the projected convergence map 
provides a high-fidelity image of the true matter distribution, allowing the dark mat- 
ter halos of individual galaxies to be viewed directly, and giving a wealth of statistical 
and morphological information about the relative distributions of mass and light. For 
instrumental designs like that planned for the Square Kilometer Array (SKA), high- 
fidelity mass imaging may be possible near the resolution limit of the core array of the 
telescope. 

Key words: large-scale structure of Universe - dark matter - gravitational lensing 
- intergalactic medium 



1 INTRODUCTION 

Since Zwicky (1933) first realised that unseen material is 
needed to explain the dynamics of galaxy clusters, many 
observations have indicated that large-scale structures are 
dominated by some form of dark matter. The now widely 
accepted cold dark matter (CDM) model provides a consis- 
tent explanation for cosmic microwave background (CMB) 
fluctuations, for type la supernova distances, for clustering 
measures from galaxy redshift surveys, for galaxy cluster 
abundances and their evolution, and for the statistics both of 
weak gravitational lensing and of Lya absorption in quasar 
spectra. The universe apparently contains about five times 
as much dark matter as ordinary baryons, providing in to- 
tal about a quarter of the closure density. According to the 
CDM picture, every galaxy has its own dark matter halo, 
which may be partially disrupted in a group or cluster to 
produce a common halo. These structures can be predicted 
in great detail by numerical simulations, but the predictions 
are yet to be convincingly verified because we are unable to 
map the dark matter distribution in enough detail to make 
a proper comparison. Weak gravitational lensing of distant 
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galaxies has allowed progress to be made, but only near the 
centres of the largest galaxy clusters is the signal-to-noise 
sufficient for true mapping. The limitations of this approach 
are clearly indicated by the recent image of a representa- 
tive field made using deep Hubble Space Telescope data by 
Massey et al. (2007b). The resolution and sensitivity of lens- 
ing maps based on galaxies are fundamentally limited by the 
finite number density and the intrinsic ellipticities of the 
sources. In this paper, we demonstrate that much higher fi- 
delity and resolution can be achieved if future observations 
allow pregalactic HI to be used as the gravitationally lensed 
source. 

The spin temperature of neutral hydrogen during and 
before the epoch of reionisation (8 ^ z ^ 300) fell out of 
thermal equilibrium with the CMB radiation, resulting in 
the absorption and emission of 21-cm radiation. There has 
been a great deal of interest in the prospect of detecting 
and mapping this radiation using radio telescopes now un- 
der construction or in planning (see Furlanetto, Oh k, Briggs 
2006, for an extensive review). This radiation provides an ex- 
cellent source for gravitational lensing studies. Structure is 
expected in the 21-cm emission down to arcsecond scales, 
and at each point on the sky there will be ~ 1000 statis- 
tically independent regions at different redshifts, and thus 



2 Hubert, Metcalf & White 



frequencies, that could in principle be observed. Gravita- 
tional lensing will coherently distort the 21-cm brightness 
temperature maps at these different frequencies. For each 
frequency, the gradient in the brightness temperature may 
be used to obtain an estimate of the lensing distortion. Since 
the intrinsic structure in the HI gas that acts as noise on 
the estimate is uncorrelated for maps at (sufficiently) dif- 
ferent frequencies, the coherent distortion of the brightness 
temperature can be measured with high accuracy if enough 
independent redshifts are observed. In this way, a map of 
the foreground matter density can be constructed (Zahn & 
Zaldarriaga 2006; Metcalf & White 2007). 

Observing the 21-cm radiation at high redshift is chal- 
lenging. Foregrounds of atmospheric, galactic, or extragalac- 
tic origin (e.g. synchrotron radiation from electrons) domi- 
nate over the 21-cm signal in the relevant frequency range 
(Furlanetto, Oh & Briggs 2006). The foregrounds are ex- 
pected to vary much less with frequency (and with position 
on the sky for galactic foregrounds) than the 21-cm signal 
from high-redshift HI structures (Zaldarriaga, Furlanetto, & 
Hernquist 2004; Zahn & Zaldarriaga 2006). Therefore, it is 
hoped that the foreground radiation and the 21-cm signal 
may be separated by modelling the foregrounds as slowly 
varying functions of frequency. 

Subtracting foregrounds from the observed radiation 
will be complicated and will contribute noise to the tempera- 
ture map. In addition to the noise from foreground residuals, 
there is irreducible noise in the mass map constructed from 
the 21-cm-lensing signal, which comes from the unknown 
intrinsic structure of the 21-cm brightness temperature dis- 
tribution. This noise cannot be reduced by increasing the 
collecting area of the telescope, by increasing the integra- 
tion time or by improving the removal of foregrounds. Met- 
calf & White (2007) showed that if the signal-to-noise in 
the brightness temperature map at each frequency is greater 
than one, then the noise in the mass map will be close to the 
irreducible value. Increasing the frequency resolution of the 
radio observations increases the number of effectively inde- 
pendent regions along the line-of-sight until the bandwidth 
becomes smaller than the radial correlation length of struc- 
ture in the brightness temperature distribution. If the band- 
width is matched to the correlation length, the irreducible 
noise is minimised. The correlation length in turn depends 
on beam size, and is smaller for smaller beams. Thus unlike 
galaxy lensing surveys, the irreducible noise decreases with 
increasing resolution for 21-cm lensing. In practise, there is 
a trade-off because smaller bandwidth means less flux, but 
this can be compensated by increasing collecting area and/or 
integration time. In this paper, we study what is achievable 
with an idealised radio telescope, so we assume that the ir- 
reducible noise level is reached, using values calculated by 
Metcalf & White (2007) as a function of beam-size and fre- 
quency. 

At least in principle, a 21-cm lensing survey will be 
much less noisy than surveys using galaxies (because of the 
larger effective number of sources) and will have a substan- 
tially stronger signal (because of the greater distance of the 
sources and the additional structure that lies in front of 
them). For a galaxy shear map, the noise increases with 
decreasing smoothing because fewer galaxies are used to es- 
timate the shear at each point of the map. The opposite is 
true for 21-cm lensing, where a smaller beam allows one to 



observe more independent sources along each line-of-sight. 
It should therefore be possible to make high-resolution im- 
ages of the matter distribution at high signal-to-noise using 
21-cm lensing, while the smallest scale over which galaxy 
lensing can map with S/N > 1 \s ^1 arcmin, even using 
an ambitious dedicated space telescope. In addition, 21-cm 
lensing will provide information about the mass distribu- 
tion at redshifts much higher than can be probed by galaxy 
lensing. The results we show below illustrate these points 
clearly. 

This paper is organised as follows. In Sec. [2] the rel- 
evant elements of lensing theory are introduced, and the 
parameters of our idealised surveys are discussed. Our lens- 
ing simulation method is described in Sec. |3] The results of 
our simulations are presented in Sec. |4] Sec. [5] contains our 
conclusions. 



2 LENSING PRELIMINARIES 

Gravitational lensing shifts the observed position of each 
point in the image of a distant source. Take the observed 
angular position on the sky to be and the position in 
the absence of lensing to be /3. The first-order distortion in 
the image is expressed by the derivatives of the mapping 
between these angles. The distortion matrix is commonly 
decomposed into the convergence n and two components of 
shear, 7, defined by 
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To lowest order and to an excellent approximation (Vale 
& White 2003), the convergence is related directly to the 
distribution of matter through 



(2a) 



AttG 



'^g{zi){T.i{e)-p{zi)5li) = Y,Ki (2b) 



with 



, ^ x-i (°° I I i\ D(z,0)D{z',z) 

g{z) = {l + z) 'J^ dz'w{z') ^^^(^,^0) • (2c) 

Here D{z\z) is the angular size distance between the 
two redshifts, and 6{x,z) is the fractional density fluctu- 
ation at redshift z and perpendicular position x. The func- 
tion E(z) = ^Qm{l + z)^ + Qa + (1 - - Qa)(1 + z)2, 
where Qm and ^^a, are the densities of matter and the cos- 
mological constant measured in units of the critical density. 
The weighting function for the source distance distribution, 
w{z), is normalised to unity. 

Equation [2b] is the multiple-lens-plane approximation, 
in which Zi is the redshift of the iih lens plane, (6) is its 
surface density, Sh is its proper thickness, and p{zi) is the 
average matter density of the universe. This approximation 
is well justified if the planes are thin compared to the range 
in redshift over which g(z) varies and Ki ^ 1 for no more 
than one lens plane. This second requirement is well justified 
for all but a very small fraction of the sky where multiple 
galaxy clusters happen to overlap in projection. 

When we consider galaxy surveys, we model the redshift 
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distribution of usable galaxies as 
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is set to obtain the median redshift Zmed appropriate for 
each specific survey (Small, Ellis & Fitchett 1994). We esti- 
mate the smoothed convergence distribution for a Gaussian 
smoothing kernel defined by 
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where denotes the angular separation between two points 
on the sky, and the 'beam diameter' A quantifies the spatial 
scale of the smoothing. For the kernel Q, the correlation 
function for the noise in the convergence map is given by 
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where rig is the number density of source galaxies on the 
sky, and cTg is the standard deviation in the magnitude of 
their ellipticities (van Waerbeke 2000). A realistic value is 
cTe = 0.3 [for example, ae = 0.32 + 0.0014(mag - 20)^ for the 
ACS camera on HST (Massey et al. 2007a)]. The proposed 
satellite SNApj^is expected to survey ^ 2% of the sky with 
an estimated galaxy density of rig lOOarcmin"^ and a 
median redshift Zmed 1.23. [For comparison, Massey et al. 
(2007b) were able to use 71 galaxies per square arcminute in 
the HST COSMOS survey.] The DUNeQ sateUite proposes 
surveying ^ 50% of the sky with rig c^i 35arcmin~^ and a 
median redshift Zmed ^ 0.9. Several proposed ground based 
surveys - LSST[3 PanSTARR^ VISTAQ " hope to reach 
source number densities comparable to DUNE. In the follow- 
ing, we adopt these two sets of parameters as our optimistic 
assessment of the parameters defining future space- and 
ground-based gdld.-xry surveys. For a Gaussian smoothing with 
A = 1 arcmin, they yield a normalization, ctn = ^/^n(^ = 0), 
of 0.012 and 0.02, respectively, for the noise correlation. A 
lower noise level 

When simulating convergence maps derived from 21-cm 
observations, we will make the approximation w{z) — 5{z — 
Zq). This is reasonable because angular size distances are 
a weak function of source redshift over the relevant range. 
The noise in the convergence map is worked out in Met- 
calf & White (2007) under the assumption that the high- 
frequency components of maps of pregalactic HI decorrelate 
with increasing redshift separation in the same way as those 
of maps of the underlying cold dark matter distribution. In 
this case, the noise is very well approximated as a Poisson 
process smoothed by the telescope beam, which we again 
model as Gaussian. This results in the same correlations as 
in Eq. ([s]) except with a different normalisation ctn. Here, we 
adopt normalisations of 0.0042 for a A = 6 arcsecond beam 
and 0.014 for a A = 1 arcminute beam. These values are rep- 
resentative for surveys that observe the 21-cm radiation at 
redshifts around zq — 12, work close to the irreducible-noise 
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limit, cover 10 MHz in frequency, and have optimal band- 
width ^ 0.05 MHz. A beam size A = Q" is very futuristic, 
since it corresponds to a densely filled array with baselines of 
order 100 km. A beam with A = 1 arcmin might be realized 
with the planned Square Kilometer Arrajj^ (SKA) (Metcalf 
& White 2007). The assumed source redshift = 12 for the 
21-cm radiation lies within the expected obvservable redshift 
range of the SKA. 



3 SIMULATIONS 

We simulated maps of the lensing convergence using the Mil- 
lennium Simulation (Springel et al. 2005) and the multiple- 
lens-plane approximation. The Millennium Simulation is 
a very large N-body simulation of cosmological structure 
formation containing 10^*^ particles in a cubic region of 
L — 500/i~^ Mpc comoving on a side. The cosmological pa- 
rameters for the simulation are: Qrn = 0.25, Ha = 0.75, and 
a Hubble constant of /i = 73 in units of 100kms~^Mpc~^. 
The initial density power spectrum is scale- invariant (spec- 
tral index n = 1) with normalisation as — 0.9. Snapshots 
of the matter distribution were stored on disk at 64 output 
times between redshift z = 127 and 2; = 0. For ^ z ^ 1, 
the snapshots are spaced at roughly 200 Myr intervals re- 
sulting in 23 snapshots in that range. Above redshift unity, 
the snapshots are spaced approximately logarithmically in 
the expansion factor. 

For each snapshot of the simulation, we project the mat- 
ter distribution in a slice of appropriate thickness onto a 
plane and place it at the snapshot's redshift along the line 
of sight. The projected matter density on these lens planes is 
represented by meshes with a spacing of 2.5/i~^ kpc comov- 
ing. In order to reduce discreteness noise while retaining the 
high resolution of the simulation, an adaptive smoothing 
kernel is used. Before projection, the mass associated with 
each particle is distributed in a spherical cloud with Gaus- 
sian density profile and rms radius equal to half the distance 
to its 64th nearest neighbour. The projected density at each 
mesh point on the lens plane is then calculated by summing 
the contributions from each particle. 

To create convergence maps over a field of 10° x 10°, we 
shoot 36000 X 36000 light rays through the series of 50 lens 
planes which extend from z = to ^ = 9. On each lens plane, 
we calculate the projected matter density at the position of 
each ray by bilinear interpolation between surrounding mesh 
points. Using Eq. [2] the convergence Hii is then calculated 
and summed for each ray to get k{0) for the assumed source 
redshift distribution [i.e. w{z) = 5{z — 12) for the 21-cm 
emission, and w[z) given by Eq. (|3| for the galaxies]. 

The convergence maps obtained by this procedure have 
a resolution of 1 arcsecond and are essentially noise-free. 
In order to simulate maps as they would be observed, we 
add independent Gaussian noise with appropriate disper- 
sion to each pixel and smooth the maps using a Gaussian 
filter representing either the radio telescope beam or the 
required smoothing in the case of galaxy lensing. This pro- 
cedure yields convergence maps with the desired resolution 
and with noise satisfying Eq. ([5|. 
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For a 10° x 10° field, the comoving volume of the back- 
ward light cone out to z = 12 is more than 32 times that of 
the Millennium Simulation. The region out to z — 2 which 
contains most of the le using structures is still 4.5 times larger 
in comoving volume than the simulation. Every simulated 
object thus contributes to the projected mass distribution 
several times. As explained in Hilbert et al. (2007), the lens 
planes were constructed by projection along a line-of-sight 
which is tilted with respect to the principal axes of the sim- 
ulation, and as a result are periodic on a rectangular cell 
of size 1.58 x 1.66/i~^Gpc^ comoving; the periodic length 
normal to the lens planes is 5.24/i~^Gpc comoving. There 
are objects that appear multiple times at the same redshift 
when z > 1.2, but the number of these cases becomes signifi- 
cant (i.e. exceeds 1/3 of all objects at a given redshift) in our 
field only for z > 2. Objects also appear multiple times at 
different redshifts. However, objects are projected on top of 
their own images only in very few (special) directions and for 
very widely separated redshifts Generally multiple copies 
of objects are almost always seen at different redshifts and 
are almost always projected onto different foregrounds and 
backgrounds. As a result, there is effectively no duplication 
of projected structure within this field, despite the fact that 
the total comoving length of the line-of-sight out to z — 12 
is more than 14 times the side of the computational box 



4 RESULTS 
4.1 Images 

Figures [l] |2] and [3] show simulated convergence maps. The 
20^ X 20^ field shown in these examples is only a small frac- 
tion (1/900) of the full 10° x 10° map that we simulated. The 
latter is too large to be displayed in sufficient detail. This 
particular field was chosen because it has a prominent mass 
concentration at upper right that is large enough to be de- 
tected in all the cases we investigate. This halo is at redshift 
z = 0.09 and has a virial mass M200 = 6.9 x 10^^/i~^ M©. 
Thus it represents a galaxy cluster similar to the Coma clus- 
ter. For this reason, the field is not typical. From our Mil- 
lennium Simulation data, we find that there is only a 3% 
probability for a random field this size to contain a halo 
more massive than 5 x 10^^/i~^ M©. 

The second largest mass concentration visible in the 
field (at lower left) is at redshift z — 1.9 and has a mass 
of 1.5 X 10^^/i~^ M0. This is also a relatively unusual event. 
From our data, we expect one halo with M200 > lO^'^/i"^ M© 
and z > 1.5 per square degree, corresponding to a 12% prob- 
ability for a 20^ x 20^ field. There are three more halos with 
masses above lO^'^/i"^ M© visible. The most prominent of 

^ Our lens-plane geometry ensures that multiple copies at the 
same or similar redshift are separated by angles equal to or not 
much smaller than the angular scale of the simulation box at that 
redshift. These multiple copies introduce artificial correlations on 
large angular scales and decrease the statistical independence of 
well-separated parts of the simulated field, but they do not affect 
correlations on smaller angular scales, nor do they alter one-point 
statistics such as the convergence probability distribution. Since 
the artificial correlations have different angular scales at different 
redshifts, multiple copies have different foregrounds and back- 
grounds in projection. 



these (left centre) is at redshift z — 0.46 and has a mass of 
1.2 X lO^^/i~^M0. On average, we expect about two such 
clusters in each 20^ x 20^ field. 

The three convergence maps in Fig. [l] are made without 
smoothing or added noise in order to illustrate the depen- 
dence on the redshift distribution of the sources. The map 
at top left gives the expected convergence distribution (at 
1" resolution) for sources at 2; ^ 12, representing the case of 
high-redshift 21-cm lensing. The map at top right uses the 
same colour table but a different source redshift distribu- 
tion, that appropriate for an optimistic space-based galaxy 
survey. The principal impression in comparing the two is 
that there is much less structure in the 'galaxy' map. This 
reflects the lowering of the overall amplitude caused by the 
smaller source distances in the galaxy case. Averaged over 
the full 10° X 10° area, the rms value cr«, of n is 0.11 in the 
21-cm map, but only 0.03 in the galaxy map (see Tab. [l] 
for a summary of the statistical properties of our simulated 
maps). The map at lower left repeats the galaxy map, but 
now the contrast is enhanced by a factor of 1 1/3 so that the 
colour range matches that of the 21-cm map. Displayed in 
this way, the two maps look similar. Nearby objects such as 
the most massive cluster appear stronger in the galaxy map, 
whereas more distant objects appear stronger in the 21-cm 
map. There are a few large structures that appear in the 
21-cm map, but are absent from the galaxy map. These are 
objects that lie beyond the redshifts assumed for the galax- 
ies. In particular, the large mass concentration at z = 1.9 is 
clearly visible in the 21-cm map, but is virtually absent in 
the galaxy map. For the reader's orientation, the lower right 
map indicates masses and redshifts for all halos in the field 
with M200 > 10^^/i~^ M©. 

The two maps in Fig.[2]illustrate expectations for 21-cm 
lensing reconstructions based on a (futuristic) radio tele- 
scope with a Gaussian beam of width A = (correspond- 
ing to baselines 100 km). The left image includes beam- 
smearing but excludes noise. A comparison with the l" 
map at the top left of Fig. ^ shows that very little detail 
is lost, and over the full 10° x 10° area the rms value of 
K is reduced by only 9% to 0.098. The right image also 
includes noise assuming the irreducible level expected for 
observations with the optimal band width for a telescope 
beam of this size (^ 0.05 MHz). This has virtually no effect 
on the image, demonstrating that such a (very large) tele- 
scope could produce high resolution mass maps with very 
high fidelity. All structures with M200 ^ 10^^/i~^ M© (these 
are indicated in Fig. [Tjl) can be clearly identified out to 
high redshift, and even many smaller halos down to masses 
M200 ^ 10^ M© are visible. The signal-to- noise at the 
scale of the beam is very high even in low density regions, so 
substantial departures from optimal conditions could be tol- 
erated without significant degradation of the resulting maps. 

Figure [3] shows maps smoothed with a Gaussian of 
width \— 1' . The colour scale is the same in all of them and 
differs from those of Fig. ^ an Fig. [2] Now, the resolution is 
similar to that obtainable with the planned Square Kilome- 
ter Array (Metcalf & White2007). The top two images are 
for sources ai z ^ 12 with no noise (left) and with noise at 
the irreducible level expected for observations at the opti- 
mal bandwidth for this beam-size (right). Again the fidelity 
of the image is high (although some differences can be seen 
in low K areas), and many of the more massive halos indi- 
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Figure 1. Convergence maps in a field of 20' x 20^ (a), (b) and (c) are all unsmoothed maps at V resolution, (a) has the distance 
weighting appropriate for HI sources at z = 12, while (b) and (c) are weighted as appropriate for a space-based galaxy survey with median 
redshift 2;med = 1.23. Both (a) and (b) follow the colour scale indicated by the bar. (c) is identical to (b) but with contrast enhanced by 
a factor of 11/3 to allow closer comparison with (a). Panel (d) identifies redshifts (colour-coded) and virial masses (size-coded) for all 
objects in the field with M200 > 10^^ M©. The eight most massive objects are labelled explicitly. 



cated in Fig.[T]l are detected individually. The middle row of 
maps are for our optimistic space-based galaxy lensing case, 
again without and with noise, while the bottom row gives the 
corresponding maps for our optimistic ground-based survey 
parameters. As at higher resolution, one is struck by how lit- 
tle structure is visible in these maps compared to the 21-cm 
case. The rms value of k over the full 10° x 10° field is smaller 
by factors of 3 and 4 in the space- and ground-based galaxy 
cases in comparison with the 21-cm case (see Table 1). 

The fidelity of the 'observed' (i.e. noisy) maps is low in 
the galaxy lensing cases. A few of the structures seen in the 
noiseless maps are still visible in their noisy counterparts, in 
particular the largest object, but many of the low- amplitude 
peaks in these maps are due to noise. In effect, only the large 
cluster at ^ = 0.09 is unambiguously detected for both sur- 
veys, while the two 10^^h~^ Mq halos at 2; = 0.46 also stand 



out above the noise in the 'space-based' map. The larger 
halo at z = 1.9 remains unseen. The fidelity of these maps 
could be improved by increasing the smoothing length, but 
this would be at the expense of losing all the individual ha- 
los. In practice, an adaptive smoothing method such as a 
maximum entropy scheme would probably be used in order 
to remove low- significance features. This would leave rather 
little structure visible in our 20^ x 20^ patch, only the highest 
peak in the 'ground-based' case. Only a few percent of fields 
this size would contain an object massive enough to be de- 
tected with high significance in a survey with these parame- 
ters. This limitation is quite evident in current ground-based 
lensing surveys (Semboloni et al. 2006; Bacon, Refregier & 
Ellis 2000; Kaiser, Wilson, & Luppino 2000; Van Waerbeke 
et al. 2000; Wittman et al. 2000). 

In order to give a better impression of the reconstruc- 
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Figure 2. 21-cm-based convergence maps for the 20^ x 20^ field shown in Fig. [TJa), but smoothed assuming a telescope beam with 
A = Whereas (a) is noise-free, noise has been added in (b) at the irreducible value for a map of this resolution. The colour scale 
indicated by the bar at right is the same as in Fig. Ilja). 



tion capabilities of a radio telescope such as SKA, Fig. [4] 
shows the map of Fig. jsj) expanded to show a 5° x 5° field. 
(Note that this is still only 1/4 of the full field we simulated.) 
Current plans for the SKA should enable this resolution to 
be reached using the dense 'core' array, but the noise level in 
the convergence map will depend on the way in which reion- 
isation proceeds. If the number density of ionised bubbles is 
large and they persist for a significant fraction of the red- 
shift range expected for SKA (7 ^ z ^ 13) then noise levels 
nearly this good can be obtained in 90 days of integration 
time. A more pessimistic scenario is that reionisation hap- 
pens very suddenly and nearly uniformly. Even if this is the 
case, and reionisation occurs near z ^ 7^ moderate fidelity 
maps at l.Sarcmin resolution should be possible and the 
same noise levels as in Fig. [3] should be attainable but on 
3 arcmin scales. In the latter case, SKA maps will be more 
noisy than Fig. [3] after 90 days of integration, although still 
of much higher fidelity than galaxy-based maps. 

The conclusion of this section is that galaxy lensing 
surveys do not provide sufficient signal-to-noise to image 
any but the most massive individual dark matter structures, 
but that a very large radio telescope could, in principle, 
provide high-resolution, high-fidelity images of the cosmic 
mass distribution in which the halos of individual massive 
galaxies and galaxy groups are visible. 

4.2 Pixel distributions 

Another useful way to represent the information in our sim- 
ulated maps is to plot the probability density function for 
the convergence, pdf(Ai:), in the different cases. For this we 
can use the full 10 x 10 degree field, rather than the smaller 
subfields discussed in section [4T] Quantitative statistics for 
all these distributions are given in Table 1. 

Figure [5] shows the results for sources at z = 12, as 
appropriate for pregalactic HI. This confirms quantitatively 
our previous conclusion that the irreducible noise has very 



little effect on the maps. Indeed, its effects are not even 
visible for a A = 6^^ beam, and they are still small for \ — l' . 
This just reflects the fact that the pdfs for the noise are 
much narrower than those for the signal, as illustrated in 
Fig. |5] Note that the narrower noise pdf is associated with 
the higher resolution, higher amplitude signal pdf. 

Figure |6] gives corresponding results for sources with 
the redshift distribution appropriate to a space-based galaxy 
lensing survey. Note that the scale has changed from Fig. [5] 
reflecting the substantially lower amplitude of fluctuations 
in n in this case. For reconstructions with a beam of width 
A = 1^, the noise expected in such a survey has a strong 
effect on pdf(Ai:). The low n tail of the observed distribution 
is practically all due to noise, and the shape of distribution 
is largely lost. Estimating the skewness or kurtosis of the 
underlying distribution would clearly require a very good 
understanding of the properties of the noise. 

The corresponding pdfs for the source redshift distribu- 
tion and noise appropriate to a ground-based galaxy survey 
are shown in Fig. [7| Here the noise destroys almost all of 
the information in the original pdf. With a large amount of 
data and with good knowledge of the systematics one can 
recover the variance accurately, but determination of higher 
moments would be extremely challenging. 

Even when the noise is high compared to the dispersion 
in At, it is still possible to measure the number density of very 
high mass objects. Figure [S] illustrates this point by plotting 
the high hi tails of the cumulative distribution functions of 
K. The noise has relatively little effect on these distributions 
for At ^ 0.1, even for the ground-based galaxy survey case. 
Such high values have a probability of around ~ 10""^ cor- 
responding to of order one object per square degree on the 
sky. For our space-based survey parameters, the noise be- 
comes unimportant for n ^ 0.05, corresponding to roughly 
100 objects per square degree. 
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Figure 3. Convergence maps smoothed at A = 1^ for the 20' x 20' field shown in Fig. ^ (a) and (b) 21-cm-based map without and with 
noise, resp.; (c) and (d) space-based galaxy lensing map without and with noise; (e) and (f) ground-based galaxy lensing map without 
and with noise. All panels use the colour scale indicated by the bar at right. 
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Figure 4. 21-cm-based convergence maps for a 5° x 5° field with a A = 1 arcmin beam and noise. The colour scale is the same as in 
Fig. [3] The field shown in Figs. ^[2] and [s] lies at the centre of this map. 



5 CONCLUSION 

The noise in a mass map constructed using gravitational 
lensing of the high-redshift HI distribution is expected to 
be much smaller than the signal. In addition, the signal-to- 
noise increases with the resolution of the map. It should thus 
be possible to make high-resolution, high-fidelity images of 
the dark matter distribution in which the dark halos of indi- 
vidual galaxies and galaxy groups are visible. For example, 
a very large future telescope with baselines ^ 100 km may 
eventually allow us to detect halos with virial masses M200 ^ 
lO^^/i-^ M© out to redshift z - 10 (Metcalf & White 2007). 



Such detailed observations will provide a very direct and 
accurate test for structure-formation models. Even with an 
SKA-like telescope, halos with M200 ^ 10^^/i~^ M© should 
be clearly detected out to high redshift. This contrasts 
strongly with mass maps constructed using gravitational 
lensing of distant galaxies, where high fidelity is only achiev- 
able for angular smoothings so large that all but the nearest 
and most massive individual objects are lost. 

Our estimates of the irreducible noise are based on a 
convergence estimator that is not necessarily optimal. It 
may therefore be possible to achieve smaller 'irreducible' 
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Figure 5. The probability density function pdf(^c) of the con- 
vergence for sources at z = 12 (black line) compared to the dis- 
tribution smoothed with a Gaussian beam with A = 6" (blue 
solid curve), and 1' (red solid curve). The dashed curves show 
the smoothed distributions with noise added at the irreducible 
level expected for observations of the pregalactic HI with opti- 
mal total bandwidth and frequency resolution. The dotted lines 
illustrate the noise distributions by themselves. 
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Figure 6. The probability density function pdf(^c) of the conver- 
gence for an optimistic space-based galaxy lensing survey (median 
redshift ^^med = 1.23). The unsmoothed, noiseless case is shown 
in black. The red curves are for Gaussian smoothing of width 
A = 1 arcmin. The dashed and solid curves are with and with- 
out noise respectively. The assumed density of source galaxies is 
100 arcmin"^. Note the difference in ^c-scale compared to Fig. [H] 

noise levels than we quote. In practice, however, it is likely 
that other sources of error will dominate the overall budget, 
for example, the error introduced by incomplete foreground 
subtraction (see Furlanetto et al. 2006, for a review). For 
most purposes, imaging the surface density does not require 
reaching the irreducible noise limit; the predicted signal is 
large enough to accommodate a noise level many times the 
irreducible value. In addition, the noise within a patch of 
area A goes down like A~^l^ while the density fluctuations 
go down roughly like so even if the noise is too large 

to map the surface density on the scale of a single beam, 
a high-fidelity map with larger smoothing can still be con- 
structed (as in the galaxy lensing case). The Square Kilome- 
ter Array in its currently proposed configuration should be 
able to map the mass distribution on arcminute scales with 
moderate to high fidelity if reionisation is not completed too 
early (Metcalf & White 2007). The optimal bandwidth for 
observing lensing is ^ 0.05 MHz while the signal-to-noise 




-0.04 -0.02 0.02 0.04 



X 

Figure 7. The probability density function pdf(^c) of the con- 
vergence for an optimistic ground-based galaxy lensing survey 
(median redshift ^^med = 0-9) • The unsmoothed, noiseless case 
is shown in black. The red curves are for Gaussian smoothing 
of width A = 1 arcmin. The dashed and solid curves are with 
and without noise respectively. The assumed source density is 
35 arcmin"^. 




Figure 8. The cumulative distribution function cdf(^t) of the 
convergence for HI at 2; '-^ 12 (red curves), for an optimistic 
space-based galaxy lensing survey (blue curves), and for an op- 
timistic ground-based galaxy lensing survey (green curves). The 
solid curves represent the distribution with a Gaussian smoothing 
of width A = 1 arcmin and with no noise. The dashed curves are 
the same but with noise included. 

for mapping the pregalactic HI at the same angular scales 
is maximal at larger band widths ^ 0.5 MHz. Lensing bene- 
fits from the stacking of many narrow redshift slices even if 
they are individually noise dominated while the temperature 
fluctuations themselves get diluted (Metcalf & White 2007). 
To reach scales of a few arcseconds as discussed here will 
require a larger telescope with dense sampling. Given the 
narrower science goals, this may be achievable with simpler 
and cheaper antennas. 

While high-resolution images of the cosmic mass distri- 
bution would be a unique product of observations of 21-cm 
lensing, they are not the only reason to carry out such stud- 
ies. If enough of the sky can be surveyed, cosmological pa- 
rameters such as the density of dark energy and its evolution 
with redshift can be measured with much higher accuracy 
by a combination of 21-cm lensing with galaxy lensing than 
they can by galaxy lensing alone or indeed by any other 
method proposed so far (Metcalf & White 2007). The base- 
line configuration of SKA may be powerful enough to achieve 
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Table 1. The expected distribution of the convergence n for 
future pregalactic HI, space- and ground-based galaxy surveys. 
For various beam sizes A and noise levels ctn, we give the standard 
deviation cr,^, the skewness S^,, the 25% quantile ^^25% 5 the 
75% quantile ^^75% of the convergence distribution. (The mean of 
1^ is zero by definition.) 



Survey 


A 




(7k 


^3 


^25% 


^75% 


21-cm, 


1'' 




0.11 


1.73 


-0.079 


0.049 


ZQ = 12 






0.098 


1.35 


-0.071 


0.046 




6'' 


0.0042 


0.098 


1.34 


-0.071 


0.046 




V 




0.058 


0.52 


-0.044 


0.031 




1' 


0.014 


0.060 


0.47 


-0.045 


0.032 


space-based, 


1" 




0.030 


3.95 


-0.020 


0.007 


Zmed — 1-23, 


1' 




0.018 


1.99 


-0.014 


0.006 


rig = lOOarcmin"^ 


1' 


0.012 


0.022 


1.09 


-0.017 


0.010 


ground-based. 


I" 




0.021 


4.61 


-0.014 


0.005 


Zmed = 0.9, 


1' 




0.014 


2.49 


-0.011 


0.004 


ria- = 35arcmin~^ 


1' 


0.02 


0.024 


0.31 


-0.018 


0.014 



much of this improvement if problems with foregrounds can 
be overcome. 
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